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Abstract 



The quark contribution to the free energy of a hot quark-gluon plasma is 
calculated to leading order in hard-thermal-loop (HTL) perturbation the- 
ory. This method selectively resums higher order corrections associated with 
plasma effects, such as screening, quasiparticles, and Landau damping. Com- 
paring to the weak-coupling expansion of QCD, the error in the one-loop HTL 
free energy is of order as, but the large a^^'^ correction from QCD plasma ef- 
fects is included exactly. 
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I. INTRODUCTION 



Experimental data from ultrarelativistic heavy-ion collisions at RHIC will soon become 
available. In order to determine if a quark-gluon plasma has been created, a careful compari- 
son of the predictions of hadronic models and QCD has to be made. It is therefore desirable 
to find a systematic way to calculate the thermodynamic properties and signatures of a 
quark-gluon plasma within QCD. Asymptotic freedom suggests that at sufficiently high 
temperatures a straightforward perturbative expansion should suffice. However, at experi- 
mentally accessible temperatures, perturbative QCD does not seem to be of any quantita- 
tive use [1-3]. The problem is evident in the free energy of the quark-gluon plasma. The 
weak-coupling expansion has been calculated through order ct'J'^ [1-3]. The successive ap- 
proximations to the free energy show no sign of converging at temperatures that are relevant 
for heavy-ion collisions. 

One possibility for improving the convergence of the perturbative predictions is to apply 
Fade approximants to the series in [Q, however, this technique can only be applied 
if several terms in the perturbation series are known. Another possibility is to use hard- 
thermal-loop (HTL) approximations within a self-consistent ^-derivable framework |p. A 
third approach is to apply HTL perturbation theory [^], which is an extension of the re- 
summation method of Braaten and Pisarski [0 into a systematic perturbative expansion. 
In two previous papers [0, we calculated the one-loop free energy of pure-glue QCD using 
HTL perturbation theory. Here we extend that calculation to include quarks, thus complet- 
ing the calculation of the free energy of a hot quark-gluon plasma to leading order in HTL 
perturbation theory. 

The paper is organized as follows. In the next section, we calculate the quark contri- 
bution to the free energy of a hot quark-gluon plasma at one loop in HTL perturbation 
theory. In section HI, we carry out the high-temperature expansion of the free energy, and 
in section IV we take the low-temperature limit. In section V, we compare the one- loop 
HTL free energy with the weak-coupling expansion of QCD. We conclude in section VI. In 
the appendix, we have collected the integrals that are required in the calculations. 

II. QUARK CONTRIBUTION TO HTL FREE ENERGY 

The one-loop HTL free energy for an S\J (Nc) gauge theory with Nf massless quarks is 

^HTL = (A^,' - 1) [{d - \)Tt + Tl\ + Nc^fT, + , (1) 

where Tt and Tl are the contributions to the free energy from transverse and longitudinal 
gluons, respectively, Tq is the contribution to the free energy from each flavor and color of 
the quarks, and AjF is a counterterm. The quark contribution is given by 

^5 = -^logdet[^-S(if)]. (2) 

The sum-integral in (|^) represents a dimensionally regularized integral over the momentum 
k and a sum over the Matsubara frequencies uJn = (2n + l)'7rT: 
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The factor of /i^ where /i is a renormahzation scale, ensures that the regularized free 
energy has the correct dimensions even for d ^ 3. The HTL quark self-energy in (EI) is 



— ^ log — 
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where k = k/A; and rriq is the thermal quark mass parameter p]. 

We proceed to calculate the quark contribution to the HTL free energy. The inverse 
quark propagator can be written as 
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where 



AoW = z.„--log— 



As{K) = k + -^\l 
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We can write the quark contribution (H) as 



Tq = -2^^\ogik 



2y^log 
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where we have separated out the free energy of an ideal gas of massless fermions. The first 
sum-integral in can be evaluated analytically. In the second sum-integral, the sum over 
Matsubara frequencies can be expressed as a contour integral: 



^ ^TT^ - [ [ duj , 



{Aq - As){Ao + A. 



k^ 
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where the contour C encloses the points uj = iuOn along the imaginary axis in Fig. 0. We 
have introduced a condensed notation for the dimensionally regularized momentum integral: 
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The integrand in (^ has logarithmic branch cuts that run from —00 to — u;+, from — a;_ 
to ci;_, from oj^ to +00, and from —k to k, where uj± are the quasiparticle dispersion relations 
that satisfy =F = 0, or 



= u;± =F 
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FIG. 1. The quark contribution to the HTL free energy can be expressed as an integral over a 
contour C that wraps around the branch cuts of log[(Ag — — k"^)]- 



UJ+ is the dispersion relation for the standard quark mode whose helicity equals its chirality. 
uj_ is the dispersion relation for the plasmino, a collective mode whose helicity is opposite to 
its chirality 0. The integrand in also has branch cuts running from —k to k due to the 
logarithms in and (|^). The contour C can be deformed to wrap around the branch cuts 
as shown in Fig. We identify the contributions from the branch cuts that end at ±uj^ 
as the quasiparticle contribution to !Fq from the quark mode. We identify the contribution 
from k < \u\ < U- a.s the quasiparticle contribution from the plasmino. The sum of these 
contributions is denoted by J^g,qp: 
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We identify the remaining contribution from \uj\ < k as the Landau-damping term and 
denote it by J-'g^Ld- 



The angle 9„ is 
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where L = log[(fc + u;)/(fc — uj)] and = k^ — uj^. The complete quark contribution to the 
free energy is the sum of the quasiparticle term ([l2|) and the Landau-damping term ([T3|) . 
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We first simplify tlie quasiparticle term. Tlie integral of cu^. is ultraviolet divergent since 
the asymptotic behavior of the dispersion relation uo^ is 



u;+{k) 



k 
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— — log 




(15) 



The integral of {lo^ — k) is convergent because the dispersion relation u)_{k) approaches the 
light-cone exponentially fast as /c — >■ oo. In order to extract the divergence analytically, we 
make a subtraction that renders the integral finite in d = 3 dimensions. The subtraction 
is then evaluated analytically using dimensional regularization. Our choice of subtraction 
integral for the quasiparticle term is 
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After subtracting this term from we can take the limit c? — > 3: 
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If we impose a momentum cutoff k < A, our subtraction integral (0) has power divergences 
proportional to A'^ and m^A^ and logarithmic divergences proportional to m^logA and 
log^ A. The quartic divergence is cancelled by the usual renormalization of the vacuum 
energy density at zero temperature. Dimensional regularization throws away the power 
divergences and replaces the logarithmic divergences by poles in o? — 3. In the limit d ^ 3, 
the individual integrals in (0) are given by ( |A.lj )- (|X7^ ) in the appendix. The result is 
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where d = 3 - 2e and = 2it'^^'^ /T{d/2). The last two integrals in (jl^) are functions of m„ 
only and must therefore be proportional to m^. Calculating the integrals numerically, their 
contributions to (p!7[) are (2.166 x 10^'^)m'^ and (—1.267 x 10^'^)m^, respectively. 

We next simplify the Landau-damping term (|13|). The temperature-independent integral 
has ultraviolet divergences from the region k ^ oo with ~ fc. We must again isolate the 
divergences by making subtractions. Our choice for the subtraction integral is 

•^iS = -2^: / /'do; 

JK JO 

Subtracting this from (|13]), we can take the limit c? — > 3: 
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If we impose ultraviolet cutoffs < A and < A on the energy and momentum, the subtrac- 
tion integral (|19|) has logarithmic divergences proportional to m^logA and m^log A. They 
cancel against the corresponding divergences in the quasiparticle subtraction integral 
The subtraction integral in (|T^ is evaluated in the limit c? — 3 using ( |A.4| )-( |Ar3| ): 
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The last integral in (^) is a function of rriq only and is therefore proportional to m^. Its 



contribution to dg) is (7.525 x 10 

Adding {^7^, (p^S]), ( pOD and (0), our final result for the quark contribution to the 
HTL free energy is 
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Since the quark contribution has no logarithmic ultraviolet divergences, the counterterm 
AjF in (|l]) is the same as in the pure-glue case 0. If we had used a momentum cutoff A 



instead of dimensional regularization, we would need a counterterm proportional to m^h? 
to cancel the quadratic divergence from the quasiparticle term (p!2[). 



III. HIGH-TEMPERATURE EXPANSION 



If the temperature is much larger than the quark mass parameter m^, the quark contri- 
bution to the free energy can be expanded in powers of rUq/T. The integral in (^ involves 
two energy and momentum scales: the "hard" scale T and the "soft" scale m^. The terms 
in the high-temperature expansion can receive contributions from both scales. Dimensional 
regularization makes it easy to separate these contributions. The soft contribution is ob- 
tained by expanding the statistical factor l/(e^'^ -|- 1) in (^) in powers of u/T. Using the 
methods in , one can show that the soft contribution to J-'q vanishes with dimensional reg- 
ularization. The hard contribution is obtained by expanding the logarithm in (^ in powers 
of rriq. The first term in this expansion is 

The integrand has single poles at a; = ±k and can be evaluated using the residue theorem. 
The momentum integrals can be evaluated analytically and in the limit ci — > 3 we obtain 
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The second term in the expansion is 
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Using the residue theorem and collapsing the contour onto the branch cuts from the loga- 
rithms, fEsD reduces to 
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The double integral can be evaluated by first integrating over k and then over uj. Expanding 
around e = and keeping terms only through e'', we obtain the finite result 
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The final result for the high-temperature expansion through order is the sum of the first 
term in (P) and the terms (0) and (^7[) : 
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IV. LOW-TEMPERATURE LIMIT 



It is useful to understand the behavior of the HTL free energy in the low-temperature 
limit where T ^ with vriq fixed. In this limit, Tq is proportional to m^. The coefficient 
could be extracted directly from the final expression for JF^, but it is simpler to compute 
it from our original expression (|^) for the quark contribution to the free energy. As T — > 0, 
the sum over the discrete Matsubara frequencies ujn = (2n + 1)ttT becomes an integral over 
the continuous Euclidean energy u: 
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After rescaling the energy u — > ku, we obtain 
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where 
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and / is the complex conjugate of /. Integrating over k, we obtain 
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Expanding around d = 3, we get 
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The integral of can be evaluated analytically and is purely imaginary. It is cancelled 
exactly by the integral of p. This cancellation is in accord with the observation that the 
quark contribution has no logarithmic ultraviolet divergences. The last integral in (^) must 
be evaluated numerically. The result is 



J^q — > (2.342 X 10^2 



(34) 



which is identical to the term in our complete expression ( ^21) for the quark contribution 
to the free energy. 



V. COMPARISON WITH WEAK-COUPLING EXPANSION 

In this section we present the numerical results for the one-loop HTL free energy (|1]) 
with Nc = ?) and Nf = 3. The quark term is given in (0). The gluon term is given in 
Ref. 0. It depends on the gluon mass parameter rrig and on a renormalization scale /X3 
associated with a logarithmically divergent integral over the three-momentum. We use the 
weak-coupling limits of the gluon and quark mass parameters ||^: 

2 27T (6 + Nf) 2 /o.X 

TUg = ^a,(/i4)T , (35) 

ml = Y^^^if^'^)'^^ ^ (36) 

where /i4 is the renormalization scale for the running coupling constant. We use a parame- 
terization of as{fii) that includes the effects of two- loop running: 



47r / 2/5ilogZ\ 



where (3o = 11 — \Nf, /5i = 51 — yA^/, and L = log (^/j^l/A^^ . For the relation between A-^ 
and the critical temperature Tc for the deconfinement phase transition, we use the result 
Tc = 1.05 calculated for Nf = 4 flavors of dynamical quarks [0. 
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The leading-order HTL free energy with Nf = 3 is shown in Fig It is scaled by the 
free energy of an ideal gas of quarks and gluons: 



ideal 



45 



21 
32^ 



Nf]T^. 



(38) 



To illustrate the sensitivity to the choices of the renormalization scales, we take their central 
values to be fi^ = 0.717 and = 27tT and we allow variations by a factor of two. 
The shaded band indicates the resulting range in predictions. The range of J^htl comes 
predominantly from variation in ^4 at the highest temperatures shown and from variations 
in /i3 at the lowest temperatures shown. With our expressions from rUg and J-'htl diverges 
either to +00 or —00 at /i4 = Aj^, depending on whether /is is greater than or less than 
our central value of 0.717 nig. The free energy of an SU{3) gauge theory with Nf massless 
quarks has been calculated in the weak-coupling expansion through order a^^^ 



45 



^0 + -^2— + 



TT 



a 



3/2 
-) + 

TT / V TT 



+ •^5 (^j + 0(a,Moga, 



(39) 



The coefficients in this expansion with /i4 = l-nT are 
•^o = l + iiV/ 
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3/2 



237.2 + 15.97A^/ - 0.413A^; + 



135 



7 



1 + \N^ \o. 



(1 + ¥i 



1/2 



799.2 + 21.96A^/ + 1.926A^J 



(40) 
(41) 

(42) 

(43) 

(44) 



The predictions from the weak-coupling expansion with = 2'kT are compared to the 
HTL free energy in Fig. |^. The expansions of the pressure truncated after orders a^, c^^J'^^ 
Qfg, and a^J'^ are shown as the dashed lines labelled 2, 3, 4, and 5. As successive terms in 
the weak-coupling expansion are added, the predictions fluctuate wildly. In addition, the 
sensitivity to the renormalization scale /i4 increases at each successive order. Of course, 
because of asymptotic freedom, the first few terms in the weak-coupling expansion will 
appear to converge at sufficiently high temperature. However, this occurs only at enormously 
high temperatures, where all the corrections to the ideal gas are tiny. For example, for 
Nf = 3, the a^/^ correction is smaller than the a<j correction only if < 37r/128. If 
we use to extrapolate to high temperature while ignoring the effects of heavier quark 
flavors, this corresponds to a temperature T > 500 T^. 

We now compare the high-temperature expansion of the HTL free energy in with the 
weak-coupling expansion. Using the values ( p5D and (^) for the thermal mass parameters, 
we flnd that the as correction is overincluded by a factor of 12(3 + A^/)/(12-|-5A^/). The a^/^ 
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FIG. 2. The free energy for QCD with Nj = 3 quarks as a function of T/T^.. The HTL free 
energy is shown as a shaded band that corresponds to varying /X3 and by a factor of two around 
their central values. The weak-coupling expansion through orders a 
as dashed lines labelled by 2, 3, 4 and 5. 



a. 



3/2 



and ot'J'^ are shown 



correction is included exactly in the leading HTL result. The overincluded a a correction 
and the large positive a^/^ correction combine with higher order corrections in the HTL free 
energy to give a negative correction that rises slowly with T as shown in Fig. |^. 

Lattice gauge theory has been used to calculate the equation of state of a quark-gluon 
plasma with Nj = 2 and Nf = 4 flavors of dynamical quarks. These calculations 



indicate that the pressure, which is the negative of the free energy, approaches that of an 
ideal gas from below. The approach to the ideal gas is more rapid than the leading order HTL 
result. For Nf = 4, it reaches 80% of the value for an ideal gas already at T = 2.5 T^. For 
higher values of T, the leading order HTL result lies significantly below the lattice results. 
This is not of great concern, because the difference can be accounted for by the next-to- 
leading order correction in HTL perturbation theory. At next-to-leading order, there are 
two-loop diagrams and one-loop diagrams with HTL counterterms. The contributions of 
order coming from the hard momentum regions of the two-loop diagrams will reproduce 
the order-as term in the conventional perturbative series (^). The contribution from the 
HTL counterterm diagram will precisely cancel the order-a^ term in the one-loop HTL free 
energy. Thus the next-to-leading order correction to J-'uTh/^ideai must approach [10(24 -|- 
7iV/)/(32 -|- 21A^/)]a;s/7r in the limit ag 0. This has the correct sign and roughly the right 
magnitude to account for the discrepancy with the lattice results. 
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VI. CONCLUSIONS 



We have completed the calculation of the free energy of a quark-gluon plasma to leading 
order in HTL perturbation theory by calculating the quark contribution. The quark term 
has a quadratic ultraviolet divergence that vanishes with dimensional regularization, but 
it has no logarithmic ultraviolet divergences. Comparing our result to the weak-coupling 
expansions for the free energy, we find that the error is of order but the large correction 
proportional to a^J"^ is included exactly. It is therefore possible that the HTL perturba- 
tive expansion for the free energy will have much better convergence properties than the 
conventional weak-coupling expansion. To verify this, it will be necessary to extend the 
calculations of the free energy to next-to-leading order in HTL perturbation theory. 
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APPENDIX A: INTEGRALS 



In this appendix, we collect the results for the integrals that are required to calculate 
the contribution from quarks to the one-loop HTL free energy. We use dimensional regu- 
larization, so that power ultraviolet divergences are set to zero and logarithmic ultraviolet 
divergences appear as poles in e. In the HTL free energy, the ultraviolet divergences are 
isolated in subtraction terms that must be expanded around e = through order e°. The 
integrals required to evaluate the subtractions in the quasiparticle terms are 
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The integrals required to evaluate the subtractions in the Landau-damping terms are 
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